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CN ' Abstract 

' We develop a new fully quantum method for determination of widths for nuclear decay by proton emission 



where multiple internal reflections of wave packet describing tunneling process inside proton-nucleus radial 
barrier are taken into account. Exact solutions for amplitudes of wave function, penetrability T and reflection 
R are found for n-step barrier (at arbitrary n) which approximates the realistic barrier. In contrast to 
semiclassical approach and two-potential approach, we establish by this method essential dependence of the 
penetrability on the starting point -Rform in the internal well where proton starts to move outside (for example, 
for T^i'^Ta the penetrability is changed up to 200 times; accuracy is |r + _R — 1| < 1.5 • 10^^^). We impose 
a new condition: in the beginning of the proton decay the proton starts to move outside from minimum 
O ' of the well. Such a condition provides minimal calculated half-life and gives stable basis for predictions. 

^ , However, the half-lives calculated by such an approach turn out to be a little closer to experimental data in 

comparison with the semiclassical half-lives. Estimated influence of the external barrier region is up to 1.5 
times for changed penetrability. 
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PACS numbers: 03.65.Xp, 23.50.+Z, 27.70.+q 
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fT^ \ and reflection, half-life 

^ ! 1 Introduction 
O , . 

^\ \ Nuclei beyond tlie proton drip line are ground-state proton emitters, i. e. nuclei unstable for emission of proton 
■ from the ground state. Associated lifetimes, ranging from 10"^ sec to few seconds, are sufficiently long to obtain 
wealth of spectroscopic information. Experimentally, a number of proton emitters has been discovered in the 
mass region A « 110, 150, and 160 (see [IHl] and references in cited papers). A new regions of proton unstable 
nuclei is supposed to be explored in close future using radioactive nuclear beams. 
^ . Initially, the parent nucleus is in quasistationary state, and the proton decay may be considered as a process 

where the proton tunnels through potential barrier. In theoretical study one can select three prevailing ap- 
proaches [5| : approach with distorted wave Born approximation (DWBA) , two-potential approach (TPA) , and 
approach for description of penetration through the barrier in terms of one-dimensional semiclassical method 
(WKBA). In systematical study these approaches are correlated between themselves, while calculation of pene- 
trability of the barrier is keystone in successful estimation of gamma widths. While the third approach studies 
such a question directly, in the first and second approaches the penetrability of the barrier is not studied and 
the width is based on correlation between wave functions in the initial and final states. However, the most ac- 
curate information on amplitudes and phases of these wave functions and correspondence between them can be 
obtained from unite picture of penetration of proton through the barrier, which is used in the WKBA approach 
(up to the second order approximation). 

The main objective of this paper is to pass from semiclassical unite description of the process of penetration 
of proton through the barrier used in the WKBA approach to its fully quantum analogue, to put a fully quantum 
grounds for determination of the penetrability in this problem. In order to provide such a formalism, we have 
improved method of multiple internal reflections (MIR, see Refs. [SHlOjl generalizing it on the radial barriers 
of arbitrary shapes. In order to realize this difficult improvement, we have restricted ourselves by consideration 
of the spherical ground-state proton emitters, while nuclear deformations are supposed to be further included 
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by standard way. This advance of the method never studied before aUows to describe dynamicaUy a process of 
penetration of the proton through the barrier of arbitrary shape in fully quantum consideration, to calculate 
penetrability and reflection without the semiclassical restrictions, to analyze abilities of the semiclassical and 
other models on such a basis. 

This paper is organized in the following way. In Sec. [21 formalism of the method of multiple internal 
reflections in description of tunneling of proton through the barrier in proton decay is presented. Here, we 
give solutions for amplitudes, define penetrability, width and half-life. In Sec. [31 results of calculations are 
confronted with experimental data and are compared with semiclassical ones. Here, using the fully quantum 
basis of the method, we study a role of the barrier shape in calculations of widths in details. In particular, 
we observe essential influence of the internal well before the barrier on the penetrability. We discuss shortly 
possible interconnections between the proposed approach and other fully quantum methods of calculation of 
widths. In Sec.[4l we summarize results. Appendixes include proof of the method MIR and alternative standard 
approach of quantum mechanics used as test for the method MIR and for the results presented. 

2 Theoretical approach 

An approach for description of one-dimensional motion of a non-relativistic particle above a barrier on the basis 
of multiple internal reflections of stationary waves relatively boundaries has been studied in number of papers 
and is known (see |llfll3j and references therein). Tunneling of the particle under the barrier was described 
successfully on the basis of multiple internal reflections of the wave packets relatively boundaries (approach 
was called as method of multiple internal reflections or method MIR, see Refs. [BHS]). In such approach it 
succeeded in connecting: 1) continuous transition of solutions for packets after each reflection, total packets 
between the above-barrier motion and the under-barrier tunneling; 2) coincidence of transmitted and reflected 
amplitudes of stationary wave function in each spatial region obtained by approach MIR with the corresponding 
amplitudes obtained by standard method of quantum mechanics; 3) all non-stationary fluxes in each step, are 
non-zero that confirms propagation of packets under the barrier (i. e. their "tunneling"). In frameworks of such 
a method, non-stationary tunneling obtained own interpretation, allowing to study this process at interesting 
time moment or space point. In calculation of phase times this method turns out to be enough simple and 
convenient [lOi . It has been adapted for scattering of the particle on nucleus and a-decay in the spherically 
symmetric approximation with the simplest radial barriers [51[71[H] and for tunneling of photons [TlllO) . However, 
further realization of the MIR approach meets with three questions. 1) Question on effectiveness. The multiple 
reflections have been proved for the motion above one rectangular barrier and for tunneling under it [TlllO l fTS ] . 
However, after addition of the second step it becomes unclear how to separate the needed reflected waves 
from all their variety in calculation of all needed amplitudes. After obtaining exact solutions of the stationary 
amplitudes for two arbitrary rectangular barriers [B]|9j , it becomes unclear how to generalize such approach for 
barriers with arbitrary complicate shape. In Ref. jl4j multiple internal reflections of the waves were studied for 
tunneling through a number of equal rectangular steps separated on equal distances. However, the amplitudes 
were presented for two such steps only, in approximation when they were separated on enough large distance, 
and these solutions in approach of multiple internal reflections were based of the amplitudes of total wave 
function obtained before by standard method (see Appendix A, eqs. (7), (18), (19) in this paper). So, we come 
to a serious unresolved problem of realization of the approach of multiple reflections in real quantum systems 
with complicated barriers, and clear algorithms of calculation of amplitudes should be constructed. 

2) Question on correctness. Whether is interference between packets formed relatively different boundaries 
appeared? Whether does this come to principally different results of the approach of multiple internal reflections 
and direct methods of quantum mechanics? Note that such interference cannot be appeared in tunneling through 
one rectangular barrier and, therefore, it could not visible in the previous papers. 

3) Question on uncertainty in radial problem. Calculations of half-lives of different types of decays based on 
the semiclassical approach are prevailing today. For example, in Ref. [15j agreement between experimental data 
of a-decay half-lives and ones calculated by theory is demonstrated in a wide region of nuclei from -'^'^^Te up to 
nucleus with Ad — 266 and Zd — 109 (see Ref. [16 for some improved approaches). In review methodology 
of calculation of half-lives for spontaneous-flssion is presented (see eqs. (21)-(24) in p. 321). Let us consider 
proton-decay of nucleus where proton penetrates from the internal region outside with its tunneling through the 
barrier. At the same boundary condition, reflected and incident waves turn out to be defined with uncertainty. 
How to determine them? The semiclassical approach gives such answer: according to theory, in construction of 
well known formula for probability we neglect completely by the second (increasing) item of the wave function 
inside tunneling region (see Ref. eq. (50.2), p. 221). In resuh, equality T'^ + R"^ ^ 1 has no any sense 
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(where T and R are coefficients of penetrabiiity and refiection). Condition of continuity for tfie wave function 
and for total flux is broken at turning point. So, we do not find reflection R. We do not suppose on possible 
interference between incident and reflected waves which can be non zero. The penetrability is determined by 
the barrier shape inside tunneling region, while internal and external parts do not take influence on it. The 
penetrability does not dependent on depth of the internal well (while the simplest rectangular well and barrier 
give another exact result). But, the semiclassical approach is so prevailing that one can suppose that it has 
enough well approximation of the penetrability estimated. It turns out that if in fully quantum approach to 
determine the penetrability through the barrier (constructed on the basis of realistic potential of interaction 
between proton and daughter nucleus) then one can obtain answer "no" . Fully quantum penetrability is a 
function of new additional independent parameters, it can achieve essential difference from semiclassical one (at 
the same boundary condition imposed on the wave function). This will be demonstrated below. 



2.1 Decay with radial barrier composed from arbitrary number of rectangular 
steps 

Let us assume that starting from some time moment before decay the nucleus could be considered as system 
composite from daughter nucleus and fragment emitted. It's decay is described by a particle with reduced mass 
m which moves in radial direction inside a radial potential with a barrier. We shall be interesting in the radial 
barrier of arbitrary shape, which has successfully been approximated by finite number N of rectangular steps: 

at i?i„in < r < ri (region 1), 
at ri < r < r2 (region 2), 




(1) 

at r^v-i <r < i?,nax (region N), 

where Vi are constants (i = l...A^). We define the first region 1 starting from point i?min, assuming that 
the fragment is formed here and then it moves outside. We shall be interesting in solutions for above barrier 
energies while the solution for tunneling could be obtained after by change i^i h. A general solution of the 
wave function (up to its normalization) has the following form: 

ij{r,e,^)^^Yi^{0,ip), (2) 



e^k.r e-*'=l^ at iJ^in < r < ri (region 1), 

a2 e'''^'' + 132 e~*'=2'', at ri < r < ra (region 2), 

(3) 

a„_ie*'="-i'' +/3Ar_ie-''="~i'', &l tn^i < r < tn-i (region N-1), 
At e*'="^ at tn^i < r < i?„iax (region N), 



where aj and /3j are unknown amplitudes. At and An are unknown amplitudes of transmission and reflection, 
YhniO, f) is spherical function, ki = ^^y2m{E — Vi) are complex wave vectors. We shall be looking for solution 
for such problem in approach of multiple internal reflections (we restrict ourselves by a case of orbital moment I = 
while its non-zero generalization changes the barrier shape which was used as arbitrary before in development 
of formalism MIR and, so, is absolutely non principal). 

According to the method of multiple internal reflections, scattering of the particle on the barrier is considered 
on the basis of wave packet consequently by steps of its propagation relatively to each boundary of the barrier 
(idea of such approach can be understood the most clearly in the problem of tunneling through the simplest 
rectangular barrier, see [TUHlIin] and Appendix A where one can flnd proof of the method and analysis of its 
properties). Each step in such consideration of propagation of the packet will be similar to on from the first 
2A'' — 1 steps, independent between themselves. From analysis of these steps recurrent relations are found for 
calculation of unknown amplitudes A*^"^ , 5'^"^ , a*-"-* and /S*^"^ for arbitrary step n, summation of these amplitudes 
are calculated. We shall be looking for the unknown amplitudes, requiring wave function and its derivative to 
be continuous at each boundary. We shall consider the coefficients Tj^, T^, and Rf, Rf, Rf as additional 
factors to amplitudes e^'^*'^^. Here, bottom index denotes number of the region, upper (top) signs "+" and "— " 
denote directions of the wave to the right or to the left, correspondingly. At the first, we calculate T^, 
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.T^_i and Rf,Rt...R 



N-l- 

rp+ ^ gi(fcj-fcj + i)rj J,- ^ ^^3 + 1 ^i{kj-kj+i)rj 

^ kj + kj + i ' ^ fcj + kj + i 

kj + fcj+i ' kj + kj+i 

Using recurrent relations: 

_ / +00 ^ \ 



(4) 



m—1 



^ m=l ' i — njj^^tij 



1 - 

and selecting as starting the following values: 



^jv-i — -^jv-i' ^1 —^11 "^1 — Ti, (6) 

we calculate successively coefficients R^_2 ■ ■ ■ Rt > ^2 ■ ■ ■ ^jv-i ^^^^ -^2^ ■ • ■ -^jv-i- finishing, we determine 
coefficients Pj-. 

the amplitudes of transmission and reflection: 

AT = f+_,, Ar^R+ (8) 
and corresponding coefficients of penetrability T and reflection R: 

Tmir — 7— I At I , Rmir = I \ . (9) 

We check the property: 

^" |^T|' + |v4K|' = f or Tmir + Rmir = 1. (10) 



ki 

which should be the test, whether the method MIR gives us proper solution for wave function. Now if energy 
of the particle is located below then height of one step with number to, then for description of transition of this 
particle through such barrier with its tunneling it shall need to use the following change: 

km^i^m- (11) 

For the potential from two rectangular steps (with different choice of their sizes) after comparison between the 
all amplitudes obtained by method of MIR and the corresponding amplitudes obtained by standard approach 
of quantum mechanics, we obtain coincidence up to first 15 digits. Increasing of number of steps up to some 
thousands keeps such accuracy and fulfillment of the property (ITU)) (see Appendix |B] where we present shortly 
the standard technique of quantum mechanics applied for the potential (|19p and all obtained amplitudes). This 
is important test which confirms reliability of the method MIR. So, we have obtained full coincidence between 
all amplitudes, calculated by method MIR and by standard approach of quantum mechanics, and that is way 
we generalize the method MIR for description of tunneling of the particle through potential, consisting from 
arbitrary number of rectangular barriers and wells of arbitrary shape. 

2.2 Width r and half-life 

We define width F of the decay of the studied quantum system by following the procedure of Gurvitz and 
Kalbermann |19j : 

V-P,Ft.T, (12) 
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where Pp is the preformation probability and F is the normahzation factor. T is the penetrability coefficient 
in propagation of the particle from the internal region outside with its tunneling through the barrier, which we 
shall calculate by approach MIR or by approach WKB. In approach WKB we define it so: 



TwKB = exp I -2 J dr\ 



(13) 



where R2 and R3 are the second and third turning points. According to |15j . the normalization factor F is 
given by simplified way by Fi or by improved way by F2 so: 





k{r') dr' - ^ 
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dr \ . (14) 



■--Ri 

The half-life r of the decay is related to the width F by well known expression: 

r = ?iln2/r. (15) 

For description of interaction between proton and the daughter nucleus we shall use the spherical symmetric 
proton-nucleus potential (at case / = 0) in Ref. [20] having the following form: 

V{r, /, Q) = vc{r) + wjv(r, Q) + vi{r), (16) 

where vc(^)j VNi^iQ) and w;(r) are Coulomb, nuclear and centrifugal components 

VRiA,Z,Q) l{l + l) / forr>r„, 

1 + exp — - — i3-^h for r<r„, 

d K 2r^ I rf^j 

Here, A and Z are the nucleon and proton numbers of the daughter nucleus, Q is the Q- value for the proton- 
decay, Vji is the strength of the nuclear component, R is radius of the daughter nucleus, is the effective 
radius of the nuclear component, d is diffuseness. All parameters are defined in Ref. [20 . Note that in this 
paper we are concentrating on the principal resolution of question to provide fully quantum basis for calculation 
of the penetrability and half-life in the problem of the proton decay, while the proton-nucleus potential can be 
used in simple form that does not take influence on the reliability of the developed methodology of multiple 
internal reflections absolutely and could be naturally included for modern more accurate models. 



3 Results 

Today, there are a lot of modern methods able to calculate half-lives, which have been studied experimentally 
well. So, we have a rich theoretical and experimental material for analysis. We shall use these nuclei: 73''Ta, 
75"'^ Re and yf^Ir. Such a choice we explain by that they have small coefficient of quadruple deformation /32 and 
at good approximation can be considered as spherical (we have ^ = 0). We shall study proton-decay on the 
basis of leaving of the particle with reduced mass from the internal region outside with its tunneling through 
the barrier. This particle is supposed to start from i?niin ^ t < Ri and move outside. Using technique of 
the and R^ coefficients in eqs. (H])-®, we calculate total amplitudes of transmission At and reflection An 
by eqs. ([8]), the penetrability coefficient Tmir by eqs. ([9]). We check the found amplitudes, coefficients Tmir 
and Rmir. comparing them with corresponding amplitudes and coefficients calculated by standard approach of 
quantum mechanics presented in Appendix [BJ We restrict ourselves by eq. (|14p for Fi and find width T by 
eq. and half- live tmir. by eq. ([TS]). We define the penetrability Twkb by eq. p^ . calculate F-width and 
half-live twkb by eqs. and ^5]) . 

3.1 Dependence of the penetrability on the starting point 

The first interesting result which we have obtained is essential dependence of penetrability on the position of the 
first region where we localize the wave incidenting on the harrier. In particular, we have analyzed how much the 
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internal boundary -Rmin takes influence on the penetrability. Taking into account that width of each interval 
is 0.01 fm, we consider point i?„iin as a starting point (with error up to 0.01 fni), from here proton begins to 
move outside and is incident on the internal part of the barrier in the first stage of the proton decay. In the 
Fig. [T] [left panel] one can see that half-live of the proton decay of vl^Ta is changed essentially at displacement 
of -Rmin- So, we establish essential dependence of the penetrability on the starting point i?start7 where the proton 
starts to move outside by approach MIR. At -Rform = 7.2127 fm this dependence allows us to achieve very close 



0,40 n 




Figure 1: Proton-decay for the 73 Ta nucleus: the left panel is for dependence of the half-life tmir on the 
starting point i?min, the central panel is for dependence on penetrability T]\jii^ on the external boundary i?max, 
the right panel is for dependence of the half-live tmir on -Rmax (here, we use i?form — 7.2127 fm where calculated 
Tmir at Rma.x — 250 fm coincides with experimental data Toxp for this nucleus). In all calculations factor F is 
the same. 

coincidence between the half-live calculated by the approach MIR and experimental data. 

3.2 Dependence of the penetrability on the external region 

The region of the barrier located between turning points R2 and R3 is main part of the potential used in 
calculation of the penetrability in the semiclassical approach (up to the second correction), while the internal 
and external parts of this potential do not take influence on it. Let us analyze whether convergence exists in 
calculations of the penetrability in the approach MIR if to increase the external boundary i?max (^^max > Ra)- 
Keeping width of each interval (step) to be the same, we shall increase i?max (through increasing number of 
intervals in the external region), starting from the external turning point i?3, and calculate the corresponding 
penetrability Tmir- In Fig. [T] [central panel] one can see how the penetrability is changed for 73^Ta with 
increasing Rmax- Dependence of the half-life tmir. on i?max is shown in the next figure [T] [right panel]. One 
can see that the method MIR gives convergent values for the penetrability and half- life at increasing of i?max- 
From such figures we find that inclusion of the external region into calculations changes the half-life up to 
1.5 times (rmin = 0.20 sec is the minimal half-life calculated at R^ < i?max < 250 fm, t^s = 0.30 sec is the 
half-life calculated at -Rmax = 250 fm, error = Tas/Tmin ~ 1-5 or 50 percents). So, error in determination of the 
penetrability in the semiclassical approach (if to take the external region into account) is expected to be the same 
as a minimum on such a basis. 

3.3 Results for the proton emitters 73'^Ta, 75^Re and 

So, the fully quantum study of the penetrability of the barrier for the proton decay give us its large dependence 
on the starting point. In order to give power of predictions of half- lives calculated by the approach MIR, we 
need to find recipe able to resolve such uncertainty in calculations of the half-lives. We shall introduce the 
following hypothesis: we shall assume that in the first stage of the proton decay proton starts to move outside 
the most probably at the coordinate of minimum of the internal well. If such a point is located in the minimum 
of the well, the half-live obtains minimal value. So, as criterion we could use minimum of half-live for the given 
potential, which has stable basis. Let us analyze which results such approach gives. We shall compare the 
half-lives calculated by approach MIR and by the semiclassical approach with experimental data. We should 
take into account that the half-lives obtained before are for the proton occupied ground state while it needs to 
take into account probability that this state is empty in the daughter nucleus. In order to obtain proper values 
for the half-lives we should divide them on the spectroscopic factor S (which we take from [5]), and then to 
compare them with experimental data. Results of such calculations are presented in Table 1. One can see that 
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Table 1: Experimental and calculated half-lives of the ground state proton emitters 73 Ta, 75 Re and 77 Ir. 
Here, is theoretical spectroscopic factor, twkb is half- life calculated by in the semiclassical approach, tmir 
is half-life calculated by in the approach MIR, twkb = twkb/ S^, tmir = tmir/ S^, Toxp is experimental 
data, i?forin is starting point in the internal well where the proton begins to move outside in the first stage 
of the proton decay, i?tp is turning point (values for S^, Tcxp are used from Table IV in Ref. [5], p. 1770; in 
calculations for each nucleus we use: i?niin = 0.11 fm, i?max — 250 fm; number of intervals in region from i?niin 
to 5 fm is 2000, in region from 5 fm to 8 fm is 500, in region from 8 fm to i?max is 5000) 



Parent nucleus 


Half-live- values, sec 


Points, fm 


Nucleus 


Q, MeV 


cth 
Op 


Twkb 


Tmir, 


Twkb 


Tmir 


^cxp 


-^form 


-R2,tp 




0.947 


0.66 


1.313- 10-1 


1.369-10-1 


1.99-10-1 


2.074- 10-1 


3.00-10-1 


3.1 


7.43 


li^Re 


1.214 


0.59 


1.5352-10-4 


1.5314- 10-* 


2.602- 10-1 


2.596-10-'^ 


3.70-10-'! 


3.32 


7.34 


ifir 


1.086 


0.51 


2.981- 10-2 


2.979-10-2 


5.85 - 10-2 


5.84- 10-2 


1.10- 10-2 


3.41 


7.46 



the calculated half-lives by MIR approach turn out to be a little closer to experimental data in comparison with 
half-lives obtained by the semiclassical approach. 

3.4 Comparison with other approaches of calculations of widths of proton decay 

Half-life of the proton decay is defined on the basis of width F which can be calculated by different approaches. 
For determination of width we shall use systematics of different approaches proposed in Ref. fS]. The proton 
emitters are narrow resonances with extremely small widths. Perturbative approach based on standard reaction 
theory could be expected to be accurate. Let us analyze two following approaches in such a direction. 

3.4.1 The distorted v^rave Born approximation method 

The resonance width can be expressed through transition amplitude, which in the distorted wave Born approx- 
imation (DWBA) is given so [5]: 

Ta+1,Z+1:A,Z = {^Ap^Ap \Vap\ (18) 

The DWBA calculations of the decay width requires knowledge of the quasistationary initial state wave function, 
^A+i, the final state wave function, "i! Ap 'ipAp, and interaction potential. The initial state wave function, "i! a+i, 
is written as a product of the daughter-nucleus wave function, and the proton wave function, The 
radial wave function of the proton ipi{r) = 5';(r)/r is found by numerically integrating the Schrodinger equation 
with one-body potential, and it should be irregular part of the Coulomb wave function, Gi (r) , in asymptotic 
region. So, such wave function is complex and it defines non-zero flux. As we use condition of continuity of 
total flux (i.e. absence of sources inside spatial region) we cannot obtain zero wave function in whole region of 
its definition, and at r = 0, in particular. 

In the final state the wave function of the decaying nuclear system can be written as a product of the intrinsic 
wave function of the proton and the daughter nucleus (an inner core). Radial part of the proton wave function 
is V'/(^) ~ Fi(r)lr, where Fi(r) is the regular Coulomb function. By other words, this wave function is real, and, 
therefore, it gives zero fiux exactly determined on the basis of the total wave function in the initial state. The 
total wave functions in the initial and final states correspond to different processes (with different total fluxes) . 
This conflrms that they, complete wave functions, do not take reflection from the barrier inside the internal 
region into account (but they are defined by different boundary conditions in the initial and final states only). 
Here, question about determination of the decay width is passed on successful determination of perturbation 
of the potential (that has another physical basis for the definition of the decay width as definition on the basis 
of the penetrability of the barrier). However, the question about separation of the total wave Junction in the 
internal region before the barrier into the incident and reflected waves remains unresolved in the D WBA method. 

Now, if we pass from real radial potential in optical model approach to complex one, then we shall introduce 
new additional independent parameter into our problem while the penetrability could be calculated for real radial 
barrier. Essential point in determination of the decay width in the DWBA method is accurate normalization of 
the wave functions in the initial and final states. It could introduce some (essential) uncertainty in calculation 
of width also while the penetrability is independent on such normalization absolutely. 



7 



One can calculate the decay width through time-reversed capture process. However, in such calculations 
shape of the barrier is approximated by inverse oscillator (or other potentials with knowing exact solutions of 
the wave function) and the penetrability for such a barrier could be calculated. It is clear that both internal 
well and external region do not take influence on results absolutely (like calculations in semiclassical approach). 
But, this is possible to resolve this problem accurately and taking whole studied shape of the potential barrier 
into account that we have demonstrated above in the fully quantum approach MIR. 

3.4.2 The two-potential method 

In the modified two-potential approach (TPA) introduced by Gurvitz and Kalbermann in |19] (details and 
examples can be found in [25], see also [SJ[1T1[23]) the decay width is defined so (see (16) in [51, and some 
details): 



where k = •\/2/x£'o/^i M is reduced mass, rs is radial coordinate of the barrier height, ipnijir) is the radial 
wave function for the first radial potential including internal well up to point tb, X/(?') is the regular radial 
wave function for the second radial potential including external region, starting from point tb and without the 
internal well and with asymptotic behavior 



Both wave functions are real and defined at different energy levels. So, in the TPA approach we do not consider 
fluxes and do not calculate penetrability. We do not study possible reflection of proton from the barrier in the 
internal well. We escape from a problem of separation of the total wave function in the internal well into the 
incident and reflected waves which takes influence on the resulting penetrability essentially (for example, for the 
simplest rectangular barrier with rectangular well such an uncorrect separation of the same exact wave function 
can give infinite penetrability that is explained by increased role of interference between incident and reflected 
waves). Success in obtaining the resulting width F is dependent on accuracy of correspondence between internal 
and external wave functions ipnij {r) and xi [i") which should be calculated from different Schrodinger equations 
with independent normalization. The correspondence between these wave functions is determined concerning 
only one boundary point (or it possible shift [13]) separating two potentials and boundary conditions at 
r = or at r — >■ oo. In contrary, the correspondence between the incident, transmitted and reflected waves 
in the MIR approach is determined concerning the barrier as the whole potential (with needed restrictions of 
the radial problem) that corresponds to fully quantum and unified consideration of penetration of the proton 
through the barrier shown in principle of non-locality of quantum mechanics. In particular, the transmitted 
wave in the MIR approach is strongly dependent on the depth of the internal well and its shape, while the 
external wave function x;(r) in the TPA approach is absolutely independent on these depth and shape (such 
a dependence can be found in the wave function (pnij (r) , but starting from the simplest WKB approach factor 
F directly includes it also). By other words, we have strong correspondence between incident, reflected and 
transmitted waves in the MIR approach and a possible week correspondence between the internal and external 
wave functions in the TPA approach. This plays the essential role in calculations of the decay widths and 
explains so large difference between the essential dependence of penetrability on the starting point in the MIR 
approach and practically full absence of such a dependence in the TPA approach. 

The simplest example demonstrated why this dependence really exists and it could be not small, can be found 
in classical tasks of quantum mechanics. Let us consider definition of the penetrability in [18] (see eq. (25.3), 



where D is the penetrability, ki and k2 are wave numbers of transmitted and incident waves, i.e. concerning 
the left asymptotic part of the potential and its asymptotic right part (see Fig. 5 in [TS!, p. 103), A is the 
transmitted amplitude of the wave function. This formula demonstrates that decreasing of the left part of 
potential increases the wave number ki (as is connected with asymptotic presentations (25.1) and (25.2) of 
waves) and, so, changes the total penetrability D. Result on the essential dependence of the penetrability of 
the starting point i?foim above has the similar sense, but has been obtained concerning the realistic barrier with 
the internal well and takes into account change of the internal amplitudes also. This contradicts with a possible 
little dependence of penetrability on the shape of the internal well in the TPA approach. So, these points seem 




(19) 



X/(0) = 0, Xi{^) ^ sin(A:r — Trl/2 + 5i) at r — cxd. 



(20) 



p. 103): 



ki 



(21) 
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to be reduction of the TPA approach, and confirm that this approach does not determine the penetrability in the 
fully quantum consideration in the problem of proton decay. At the same time, comparison of results obtained 
by such approach and results obtained by principally other fully quantum developments sometimes leads to 
some confusion as the TPA approach has been called as the fully quantum. So, approaches for determination 
of the decay widths on the basis of penetrability are physically motivated, could be more accurate and have 
perspective for research. 

4 Conclusions 

The new fully quantum method (called as the method of multiple internal reflections, or MIR) for calculation 
of widths for the decay of the nucleus by emission of proton in the spherically symmetric approximation and 
the realistic radial barrier of arbitrary shape is presented. Note the following: 

• Solutions for amplitudes of wave function (described motion of the proton from the internal region outside 
with its tunneling through the barrier), penetrability T and reflection R are found by the method MIR 
for n-step radial barrier at arbitrary n. These solutions are exactly solvable and have been obtained in the 
fully quantum approach for the first time. At limit n — > oo these solutions could be considered as exact 
ones for the realistic proton-nucleus potential with needed arbitrary barrier and internal hole. Estimated 
error of the achieved results is \T + R — 1\ < 1.5 • 10~^^. 

• In contrast to the semiclassical approach and the TPA approach, the approach MIR gives essential depen- 
dence of the penetrability on the starting point inform inside the internal well where proton starts to move 
outside in the beginning of the proton decay. For example, the penetrability of the barrier calculated by 
MIR approach for yl^Ta is changed up to 200 times in dependence on position of i?form (see Fig. [1] the 
left panel) . The amplitudes calculated by MIR approach we compared with the corresponding amplitudes 
obtained (for the same potential) by independent standard stationary approach of quantum mechanics 
presented in Appendix [B] and we obtained coincidence up to first 15 digits for all considered amplitudes. 
This important test confirms that presence of the essential dependence of the penetrability of the starting 
point Rioim is result independent on the fully quantum method applied. Such a result could be connected 
with a possibility to introduce initial condition which could be imposed on proton decay in the fully quan- 
tum consideration. Comparison with the WKB and TPA approaches shows that such approaches have 
no such a perspective (having physical sense and opening a possibility to obtain a new information about 
the proton decay), which fully quantum study of the penetrability gives. 

• In order to resolve uncertainty in calculations of the half-lives caused by the dependence of the penetrability 
on i?form, we have introduced the hypothesis: in the first stage of the proton decay the proton starts to 
move outside at the coordinate of minimum of the internal well. Such condition provides minimal value for 
the calculated half-life and gives stable basis for predictions in the MIR approach. However, the half-lives 
calculated by the MIR approach turn out to be a little closer to experimental data in comparison with 
the half-lives obtained by the semiclassical approach (see Tabl. 1). 

• Taking the external region of the potential after the barrier into account, half-live calculated by the MIR 
approach is changed up to 1.5 times (see Fig. [1] the right panel). 

A main advance of the MIR method developed in this paper is not a new attempt to describe experimental data 
of half-lives more accurately than other approaches do this, but rather this method seems to be the first tools for 
estimation of the penetrability of any desirable barrier of the proton decay in the fully quantum consideration. 

A Tunneling of packet through one-dimensional rectangular step 

Main ideas and formalism of the multiple internal reflections can be the most clearly analyzed in the simplest 
problem of tunneling of the particle through one-dimensional rectangular barrier in whole axis [6Hl0] . Let us 
consider a problem of tunneling of a particle in a positive z-direction through an one-dimensional rectangular 
potential barrier (see Fig. [2]). Let us label a region I for a; < 0, a region II for < a: < a and a region III for 
x> a, accordingly. We shall study an evolution of its tunneling through the barrier. In standard approach, with 
energy less than the barrier height the tunneling evolution of the particle is described using a non-stationary 
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Figure 2: Tunneling of tlie particle through one-dimensional rectangular barrier 



propagation of WP 



where stationary WF is: 



-f-oo 



il}{x,t)^ / g{E- E)ip{k,x)e~'^'/^dE, 



ip{x) 



+ Aiie-"''', for a; < 0, 



ae 



f3e 



A ^tkx 



for < a; < a, 
for a; > a 



(22) 



(23) 



and k = j-\/2mE, ^ = j-y^2m{Vi — E), E and m are the total energy and mass of the particle, accordingly. 
The weight amplitude g{E — E) can be written in the standard gaussian form and satisfies to a requirement 
of the normalization / \g{E — E)\'^dE = 1, value E is an average energy of the particle. One can calculate 
coefficients At, Ah, a and /3 analytically, using a requirements of a continuity of WF ip{x) and its derivative 
on each boundary of the barrier. Substituting in eq. ()22p instead of Lp{k,x) the incident ipinc(k,x), transmitted 
iftrik, x) or reflected part of WF iprefik, x), defined by eq. ((23|). we receive the incident, transmitted or reflected 
WP, accordingly. 

We assume, that a time, for which the WP tunnels through the barrier, is enough small. So, the time 
necessary for a tunneling of an a-particle through a barrier of decay in a-decay of a nucleus, is about 10^^^ 
seconds. We consider, that one can neglect a spreading of the WP for this time. And a breadth of the WP 
appears essentially more narrow on a comparison with a barrier breadth. Considering only sub-barrier processes, 
we exclude a component of waves for above-barrier energies, having included the additional transformation 



g{E-E)^giE~E)eiVi-E), 
where ^-function satisfies to the requirement 



(24) 



0{V) = 



0, for 77 < 0; 

1, for 77 > 0. 



The method of multiple internal reflections considers the propagation process of the WP describing a motion 
of the particle, sequentially on steps of its penetration in relation to each boundary of the barrier |11H13) . Using 
this method, we find expressions for the transmitted and reflected WP in relation to the barrier. At the first 
step we consider the WP in the region I, which is incident upon the first (initial) boundary of the barrier. 
Let's assume, that this package transforms into the WP, transmitted through this boundary and tunneling 
further in the region II, and into the WP, refiected from the boundary and propagating back in the region I. 
Thus we consider, that the WP, tunneling in the region II, is not reached the second (final) boundary of the 
barrier because of a terminating velocity of its propagation, and consequently at this step we consider only two 
regions I and II. Because of physical reasons to construct an expression for this packet, we consider, that its 
amplitude should decrease in a positive x-direction. We use only one item l3 exp{~^x) in eq. (j23p . throwing the 
second increasing item a exp(^a;) (in an opposite case we break a requirement of a finiteness of the WF for an 
indefinitely wide barrier). In result, in the region II we obtain: 



+00 



^Plix,t)= / g{E-E)e{Vi~E)l3"e~^''-'^*/^dE,ioY0 



< X < a. 



(25) 
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Thus the WF in the barrier region constructed by such way, is an analytic continuation of a relevant expression 
for the WF, corresponding to a similar problem with above-barrier energies, where as a stationary expression 
we select the wave exp(ifc2a;), propagated to the right. 

Let's consider the first step further. One can write expressions for the incident and the reflected WP in 
relation to the first boundary as follows 

+00 _ 

^i„^{x,t) = / g{E - E)0{Vi - E)e''''=-'^^/''dE, for a; < 0, 

" (26) 



/ 





iPl^j:{x,t) = / g{E - E)9{Vi - E)A%e-''''=-'^*/^dE, for a; < 0. 



A sum of these expressions represents the complete WF in the region I, which is dependent on a time. Let's 
require, that this WF and its derivative continuously transform into the WF ([25]) and its derivative at point 
X — (we assume, that the weight amplitude g{E — E) differs weakly at transmitting and reflecting of the 
WP in relation to the barrier boundaries). In result, we obtain two equations, in which one can pass from the 
time-dependent WP to the corresponding stationary WF and obtain the unknown coefficients /3'^ and A'^. 

At the second step we consider the WP, tunneling in the region II and incident upon the second boundary 
of the barrier at point x — a. It transforms into the WP, transmitted through this boundary and propagated 
in the region III, and into the WP, reflected from the boundary and tunneled back in the region II. For a 
determination of these packets one can use eq. (|22l) with account eq. ([M)) . where as the stationary WF we use: 

(pf^^{k,x) = iplr{k,x) = /3"e~'^^, for < a; < a, 

<fl{k,x) = ^Og^fea;^ forx>a, (27) 

ipl^j{k,x) — a^e^^, for < X < a. 

Here, for forming an expression for the WP reflected from the boundary, we select an increasing part of the 
stationary solution a° exp(^a::) only. Imposing a condition of continuity on the time-dependent WF and its 
derivative at point x — a, we obtain 2 new equations, from which we find the unknowns coefficients A'^ and a". 

At the third step the WP, tunneling in the region II, is incident upon the first boundary of the barrier. Then 
it transforms into the WP, transmitted through this boundary and propagated further in the region I, and into 
the WP, refiected from boundary and tunneled back in the region II. For a determination of these packets one 
can use eq. ([22l) with account eq. (f24|) . where as the stationary WF we use: 



^Lcik,x) = ipl^f{k,x), forO<a;<a, 

iptik,x) = A)je-*'=^, fora;<0, (28) 
(pl^j{k,x) = /3^e~'^'', for < a; < a. 

Using a conditions of continuity for the time-dependent WF and its derivative at point x = 0, we obtain the 
unknowns coefficients A]^ and /3^. 

Analyzing further possible processes of the transmission (and the refiection) of the WP through the bound- 
aries of the barrier, we come to a deduction, that any of following steps can be reduced to one of 2 considered 
above. For the unknown coefficients a", and A^, used in expressions for the WP, forming in result of 

some internal refiections from the boundaries, one can obtain the recurrence relations: 



fc + «| + k + k ^29) 

k + i^' ^^-^ zC-l-fc ' " ~ ^^ + k■ 

Considering the propagation of the WP by such way, we obtain expressions for the WF on each region which 
can be written through series of multiple WP. Using eq. (|22|) with account eq. (p4| . we determine resultant 
expressions for the incident, transmitted and reflected WP in relation to the barrier, where one can need to use 
following expressions for the stationary WF: 



(Pinc{k,x) = e'''^, for a; < 0, 

(ftr{k,x) ^ ^ A^e^^^ , for x > a, 

+00 

(Pref{k,x) = J2 ^Re'"'""', forx<0. 

n=0 



(30) 
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Now we consider the WP formed in result of sequential n reflections from the boundaries of the barrier and 
incident upon one of these boundaries at point a; = (i = 1) or at point x — a [i = 2). In result, this WP 
transforms into the WP ipl^{x,t), transmitted through boundary with number i, and into the WP 'il'l,^j:{x,t), 
reflected from this boundary. For an independent on x parts of the stationary WF one can write: 



- T+ 



fine 



- U+ 



exp(— ^x) ^ exgiikx) ' exp(— z/cx) 

^Ir ^ J.+ "fief ^ ^+ 

exp{ikx) ^ exp(— ^x) ' exp(^a::) ^ exp(— ^x) ' 



e:xjp{ikx) ' 

rinc 



(31) 



exjp(—ikx) 



fine 



fref 



exp(^a;) ' exp(-^a;) 



fine 

exp(^a;) ' 



where the sign "+" (or "-") corresponds to the WP, tunneling (or propagating) in a positive (or negative) 
x-direction and incident upon the boundary with number i. Using and Rf, one can precisely describe an 
arbitrary WP which has formed in result of n-multiple reflections, if to know a "path" of its propagation along 
the barrier. Using the recurrence relations eq. (|29|) . the coefficients and Rf can be obtained. 



^1 = 



An 

^2 - 



Tf = 



an+l 



(32) 



Using the recurrence relations, one can find series of coefficients a", /S", and A^. However, these series 
can be calculated easier, using coefficients and Rf. Analyzing all possible "paths" of the WP propagations 
along the barrier, we receive: 



+ 00 
+ 00 
+ 00 

E a" 

+ 00 

E r 

n=0 



-\-oo 



i4:k£_e 



Fsub 



^ T/Tf ( 1+ Y^iRtRiT] = 

\ n=l 

= Rt+ T+RtT,- (l + EiRtR^r) = ^ 

\ n=l J suh 



+ E {RtRx T ) = 

/?"(!+ E(fi2+^r; 



2k{ii~k)e 

-^sub 

2k{i£_ + k) 

^ sub 



(33) 



where 



D± 

k"^ 



(fc2 - ^2)^ 

1 ± e"2Ca 



2ik^D^ 



k^+e 



2mVi 



(34) 



All series E Z^". E and 

with the corresponding coefficients a, /3, At and Aji 
Using the following substitution 

k2 



, obtained using the method of multiple internal reflections, coincide 
of the eq. (f23l) . calculated by a stationary methods [18] . 



(35) 



where k2 — j- \j2m{E — Vi) is a wave number for a case of above-barrier energies, expression for the coefficients 
a", and for each step, expressions for the WF for each step, the total eqs. ([33]) and (p4l) transform into 

the corresponding expressions for a problem of the particle propagation above this barrier. At the transformation 
of the WP and the time-dependent WF one can need to change a sign of argument at 6'-function. Besides the 
following property is fulfilled: 



+ 00 



n=0 



+ 00 



ri=0 



1. 



(36) 



B Direct method 



We shall add shortly solution for amplitudes of the wave function obtained by standard technique of quantum 
mechanics which could be obtained if to use only condition of continuity of the wave function and its derivative 
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at each boundary, but on the whole region of the studied potentiaL At first, we find functions /2 and 52 (from 
the first boundary): 

/2 = ^e2''=-S g, = ^e'('=+'^-)-. (37) 



Then, using the fohowing recurrent relations: 

^^+^ (fc,+i+fc,)e2''=.-^+/,(fc,+i-%) ' ^ ^ 



we calculate next functions /a, /4, ■ ■ ■ fn, and by such a formula: 



2kj e'(''j+i+'^j)^J 
(fc,+i + fc^pTfe-^TT/T^j+i - kj) 



the functions 53, (74, (75 . . From /„ and g„ we find amplitudes q;„, /3„ and amplitude of transmission At'. 



0, At = an = -^. (40) 



Now using the recurrent relations: 



= ^^^^ J ^^^i^ (41) 



and such a formula: 

ft +5,, (42) 

we consistently calculate the amplitudes a„_i, /3„-i, 0:71-2, /?n-2 ■ • ■02, /32- At finishing, we find amplitude of 
reflection Aj^: 

An = a2 e'^'^+'^^^'^i + e'^'"'''^^''' - e'^'''''' . (43) 

As test we use condition: 

^|^tP + |AhP = 1. (44) 
fci 

Studying the problem of proton decay, we used such a techniques for check the amplitudes obtained previously 
by the MIR approach and obtained coincidence up to first 15 digits for all considered amplitudes. In particular, 
we reconstruct completely the pictures of the probability presented in Fig. [T] (a) and (b) , but using standard 
technique above. So, result on the large dependence of the penetrability of the position of the starting point 
Rfovm in such figures is independent on the used method. 
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